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Quantum turbulence that exhibits vortex creation, annihilation and interactions is demonstrated 
as an exact solution of the time-dependent, free-particle Schroedinger equation evolved from a 
smooth random-phased initial condition. Rela:xed quantum turbulence in 2D and 3D exhibits uni- 
versal scaling in the steady-state energy spectrum as k^^ in small scales. Due to the lack to 
dissipation, no evidence of Kolmogorov-type forward energy cascade in 3D or inverse energy cas- 
cade in 2D is found, but the rotational and potential flow components do approach equi-partition 
in the scaling regime. In addition, the 3D vortex line-line correlation exhibits universal behavior, 
scaled as Ar~^, where Ar is the separation between any two vortex line elements, in fully developed 
turbulence. We also show that quantum vortex is not frozen to the matter, nor is the vortex motion 
induced by other vortices via Biot-Savart's law. Thus, quantum vortex is actually a nonlinear wave, 
propagating at a speed very different from a classical vortex. 

PACS numbers: 67.10.Fj, 47.27.-i, 47.32.-y, 47.50.-d, 05.10.-a, 47.ll.-j 



I. INTRODUCTION 

Vortices are probably the most common features in 
classical fluid turbulence, regardless of compressible or in- 
compressible flows. Likewise, quantized vortex lines and 
tangles not only are frequently observed but also play 
vital roles in quantum turbulence [1,2,3]. Quantized vor- 
tices also exist in type II superconductors, representing 
defects in superconductivity within which resistivity is fi- 
nite [4] . Recent interests in quantum vortices have been 
motivated by the intriguing dynamics in Bose-Einstein 
Condensation [5.6,7]. Particularly interesting is the exis- 
tence of Kelvin waves on vortex lines [8,9]. Kelvin waves 
are helical waves and believed to be responsible for cou- 
pling quantum turbulence to dissipation. 

While two dimensional quantum vortices can be cre- 
ated or annihilated in pair, three dimensional vortex lines 
can change their topology through reconnection. These 
intriguing dynamical behaviors have previously been re- 
ported with a simple model [10]. However, it is not clear 
how quantum turbulence may arise from abundant quan- 
tum vortex excitations, with the dynamics of each vortex 
obeying what is described by that work. 

Turbulence can occur when an abundance of free en- 
ergy is available. External forcing is one option to pro- 
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vide the free energy, yielding the driven quantum tur- 
bulence. A configuration unstable to large scale fiuctua- 
tions, such as that due to gravity [11], is an alternative 
to make the free energy available, creating the relaxed 
turbulence. Driven turbulence requires an energy sink or 
dissipation so as to achieve detailed balance in a steady 
state. Relaxed turbulence on the other hand does not 
necessarily involve in dissipation. But when it does, tur- 
bulence will only be a transient phenomenon and will 
eventually decay away. On the other hand, dissipation- 
less relaxed turbulence that obeys Hamiltonian dynamics 
with constant energy is a rather uncommon phenomenon 
in a macroscopic world, since when all degrees of free- 
dom are excited the smallest scale in real systems can 
often couple to dissipation. Moreover in most theoretical 
models of turbulence, integration of dynamical equations 
must often involves numerical dissipation so as to war- 
rant numerical stability, thereby rendering investigation 
of dissipationless turbulence difficult. 

Dissipationless dynamics can nevertheless exist in iso- 
lated microscopic systems. How can a microscopic dy- 
namical system evolve to a turbulence state? A turbu- 
lent system is normally conceived to be chaotic in detail 
but universal (insensitive to initial conditions) on aver- 
age. It is unquestionable that a many-body system with 
nonlinear interactions can be chaotic. However, it is un- 
clear how a universal state can dynamically arise from a 
Hamiltonian system. One of the goals of this work is to 
demonstrate such a universal non-equilibrium state in a 
Hamiltonian system. 

Another important feature which may also exist in tur- 
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bulence is the coherent structure [12]. Coherent struc- 
tures are comparatively longer-hved than other more 
common random flows in turbulence, and coherent struc- 
tures may even possess correlation among different length 
scales of great disparity. The latter property suggests 
that local scale coupling, which is the basis of the Kol- 
mogorov turbulence cascade scenario, does not necessar- 
ily hold, and can lead to intermittency and deviation 
from the universal turbulence spectrum derived from the 
mean-field theory [13]. Quantum vortices may be re- 
garded as such coherent objects existing in quantum tur- 
bulence, and we will show that their existence can alter 
the turbulence small-scale energy spectrum. 

Inspired by the concept developed for classical fluid 
turbulence, we are motivated to examine quantum tur- 
bulence in the framework of vortex creation and anni- 
hilation. Due to the strong nonlinearity in turbulence 
we shall in principle investigate this problem numeri- 
cally. However, numerical simulations can often yield 
uncontrollable numerical dissipation that depends on the 
adopted numerical schemes. Because the quantum dy- 
namics obeys certain exact symmetries, such as the uni- 
tary condition, which can be difficult to be strictly rein- 
forced in numerical simulations, we therefore study a sim- 
plest, soluble but non-trivial case in this work — a system 
described by the Schrocdinger equation of free particle 
with random-phased initial conditions. In the context of 
superfluid, our system corresponds to the regime where 
the vortex healing length diverges. 

This system in wavefunction formulation is linear and 
exactly soluble. However, in the fluid formulation, this 
linear quantum system can be strongly nonlinear. When 
compared with the nonlinear Schroedinger equation sys- 
tem, the difference is simply the absence of polytropic 
ffuid pressure in the linear Schroedinger case. The poly- 
tropic fluid pressure is responsible for sound emission out- 
side vortex cores and important only on large scale. Due 
to the existence of the polytropic fluid pressure, in ad- 
dition to the forward cascade of kinetic energy the oc- 
cupation number per mode can also exhibit an inverse 
cascade in favor of large-scale excitation. (In the linear 
Schroedinger case, the occupation number per mode is 
conserved and the inverse cascade is prohibited.) But 
important features of quantum turbulence, such as vor- 
tices, are of small scale. It is therefore anticipated that 
the linear Schroedinger equation should be able to cap- 
ture these small-scale features. 

This reasoning is strengthened by the following obser- 
vation of duality for the nonlinear Schroedinger system. 
At the very location of vortex where the fluid formu- 
lation becomes strongly nonlinear and the density van- 
ishes, the system becomes linear in the wave formula- 
tion. The same argument also applies to Kelvin waves 
which are located at the low-density vortex cores. By 
contrast, the strongly nonlinear regime of a nonlinear 
Schroedinger system in the wave formulation can be ap- 
proximated in the fluid formulation by a linear equation 
for sound waves. The duality conjugation of strong cou- 



pling regimes between these two different formulations 
leads us to believe that our linear wave model is essen- 
tial to understand quantum vortices and has generic rel- 
evance to quantum turbulence of the full problem de- 
scribed by a nonlinear Schroedinger equation. 

This paper is organized as follows. Sec. (II) presents 
the fluid formulation, an alternative to the conventional 
wave formulation, and the connections of the two in terms 
of vortices. We give the numerical results of quantum 
turbulence prior to full excitation of vortices in Sec. (HI). 
Vortices in fully developed turbulence are addressed in 
Sec. (IV). Our conclusion is given in Scc.(V). Throughout 
this work, we let fi. = 1 and mass m = 1. 



II. FLUID FORMULATION 

To show the nonlinearity in the fluid variables of the 
free-particle Schroedinger equation, we cast the wave 
equation into a set of fluid equations. Let the com- 
plex wavefunction be represented by two real functions, 
^ = fe^^, with / representing the real amplitude and S 
the real phase. We further let the fluid density be p = p 
and the fluid velocity v = Im{-^*V^)/p = \/S. The 
Schroedinger equation then becomes 



at 



+ V • ifWS) = 



and 



dVS (V5)^ 
dt ^ 2 
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1=0. 
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The first equation represents the conservation of mass. 

We also replace the original V{dS/dt) by d{'VS)/dt in 
the second equation, representing the Euler equation for 
fluid momentum. 

On the other hand, when S becomes a multi-valued 
function, vortices appear and rotational flows set in. The 
rotational flow is induced by vortices via Biot-Savart 
law. The multi-valued S can appear when p locally 
vanishes and the quantum potential V^//2/ diverges. 
Therefore vortices arc always located at regions where 
the wavefunction vanishes [10,14]. Since the wavefunc- 
tion is complex, this condition demands both real part 
R and imaginary part / to vanish simultaneously, i.e., 
* = (x - xo) • VR + i(x - xo) • V/. Normally VR and 
V/ are not aligned at x = xq, and the intersection of the 
two non-parallel surfaces R = and 1 = forms a line 
in 3D, the vortex line. In 2D, the intersection is a point 
vortex. 

The local density near the vortex can be expressed 
as p = [(x — Xq) • V7?]^ + [(x — Xo) • VJ]^, which van- 
ishes quadratically. Upon diagonalization it becomes 
p = a^{x' - x'^ityf -h \?{y' - y'oit))'^ where a and b are 
constants. Without loss of generality, we can therefore 
let ~ a{x — xo(t)) + ib(y — yo(t)) in an instantaneous 
diagonal frame where the primes in the above expression 
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is dropped for convenience. Prom this local expression of 

wavefunction, we find tfiat 



* = A/a2(x-a;o)2 + 62(y_yo) 



2j<t> 



(3) 



where S = <j> which is the space rotation angle around the 
vortex line. Thus, isolated vortex is generally quantized 
with a unit angular momentum. (Vortex quantized with 
an integer multiple of unit angular momentum is also 
allowed, but it is much less common to be excited in 
dynamical evolution.) Near the vortex line the leading- 
order velocity 



V5: 



abr 



(4) 



where r = x — Xq. The velocity diverges at the vortex 
and the singular vorticity can be defined as 



V X VS' = ^ lim VS-d\ = 5{r)z, (5) 



where the subscript c represents a circular contour along 
4> and Tc is the radius of the circle c. This definition of curl 
operation makes the singular vorticity circular symmet- 
ric. In fact when the density is not circularly symmetric, 
the leading order velocity also contains compression V^S", 
which also diverges. 



V^5 = 



a6(6^ - g^) sin(2</.) 



One can further derive that 



a 



(6) 
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to the leading order near the vortex, and we can 
uniquely separate the rotational component S'^ from 
the potential component S^, with = (p and = 
tan~^[(b/a) tan(0)] — </> at the vortex. 

In addition, vortex can move at a finite velocity per- 
pendicular to the vortex line. Since the vortex motion 
is locally two-dimensional, we shall consider the 2D case 
in the following local analysis. The vortex motion can 
be determined from the density evolution equation, since 
the vortex is tied to p = 0. Substituting the local expres- 
sion of v and p into Eq.(4), we find that to the leading 
order V • (pv) = abV ■ (rcj)) = 0. Let p = po + 5p and 
V = vo -h <5v, where po and vq are the leading order den- 
sity and velocity given above. We then obtain 

V ■ {po6v + Spvo) = (r • Vi?)V^/ - (r • V/)V^i? (8) 

to the next order. Since this equation is a linear equation 
in X and y, we can regard it as the convection of the den- 
sity null with a finite velocity w; that is, w • Vpo where 
w = dxo{t)/dt, describing the instantaneous velocity of 
a vortex. A straightforward algebra shows that 



2V^' X V** • z ' 



(9) 



where the gradient V is two-dimensional perpendicular 
to the local vortex line. Note that the propagation ve- 
locity is primarily determined by the Laplacian of wave 
function. This expression is the same as that obtained in 
[10] with a somewhat different approach. 

With this expression for w, we may inquire whether w 
is the same as, or can be reduced from, the material ve- 
locity Sv. This question is important since if it does the 
vortex will be frozen to the quantum fluid and gets ad- 
vected as the matter does, much like the classical vortex 
in an isentropic classical fluid. If it does not, the quan- 
tum vortex is then a wave pattern, propagating with its 
own velocity. As far as we are aware of, this important 
issue has not been addressed in the literature. 

Since the vorticity w(s V x v) is a Dirac (5-function, 
which has circular symmetry, the zeroth-order velocity, 
being along the circular direction, can not displace the 
vortex. That is, vq • Vw = 0. It is the next-order ve- 
locity Sv that can provide a non-vanishing contribution 
to vorticity advection. Since the leading-order velocity 
is singular, the next-order velocity can only be obtained 
by regularization. Regularization of velocity can be per- 
formed by averaging the continuity equation, Eq.(l), over 
an infinitesimally small disk around the vortex: 



dt 



+ V • {pVS) = 0. 



(10) 



The continuity equation is used to defined the advection 
velocity because it ensures that the averaged material the 
vortex is advected by an effective material velocity. Such 
an average is equivalent to a local angular average, and 
the regularized velocity 6v can be shown to be 



WpV5)) 



/m[V^* • (VV - 

2|WP ' 



(11) 



after a straightforward algebra. Here the identity ma- 
trix I is two dimensional and so is the gradient V. The 
material velocity 5v near the vortex is proportional to 
the quadrupole expansion of the wave function and ap- 
parently different from the vortex velocity w which is 
proportional to the Laplacian of the wave function. It 
evidences that the vortex is not frozen to the matter and 
must therefore be a wave pattern. 

Finally, we may want to examine whether the vortex 
motion arises from the induction velocity generated by 
other vortices via the Biot-Sarvart law. Note that the 
induction velocity is incompressible, whereas the mate- 
rial velocity described by Eq.(ll) is composed of both 
compressible and incompressible components and there- 
fore the induction velocity is different from the material 
velocity. This question is important as it has often been 
assumed so in the literature. We find that the answer is 
again negative, as can be illustrated from the following 
simple case. 

Consider a two-dimensional free-particle wave func- 
tion in cylindrical coordinate {R,(j)), ^ = cq — 
Ji(fci?)e*('^~'^ where cq is a positive constant, 
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the energy eigenvalue and Ji the first-order Bessel func- 
tion. The density jvl/p is stationary in a counter- 
clockwise rotating frame with an angular velocity Rk^ /2 
around R = Q. This system contains only two vor- 
tices of opposite signs near the first peak of Ji when 
Ji{kRi) > Co and Ji(kR2) < cq, where Ri and R2 arc 
the first and second peaks of Ji. Since the two density 
nulls are also stationary in the same rotating frame, the 
two vortices travel together counter-clockwise with the 
same angular velocity. If the vortex motion were to be 
due to the induction of the other vortex, the two should 
have traveled with a speed depending on their separation. 
One may change the vortex separation by adjusting cq. 
However the angular velocity of the vortex pair remains 
the same under such a change, meaning that the vortex 
motion is not caused by the mutual induction. 

The nature of a quantum vortex is therefore very dif- 
ferent from a classical vortex of an isentropic classical 
fluid, for which the fluid dynamics is described by the 
same mass contimiity equation as Eq.(l) and a similar 
momentum equation as Eq.(2). The classical vortex is 
frozen to the matter. However, the quantum vortex is a 
nonlinear wave, whose propagation is controlled by the 
local Laplacian of wave function. As a consequence, the 
dynamics of quantum vortex, or that of the rotational 
velocity component, cannot be described by the fluid for- 
mulation given by Eqs.(l) and (2) alone for a potential 
flow. When vortices appear in a quantum fluid, an ad- 
ditional equation of motion for vortex is needed in the 
fluid formulation to evolve the whole system. 

To trace the origin of this issue, let us return to the 
momentum equation of the fluid formulation, Eq.(2). We 
shall prove from the vortex motion that the time deriva- 
tive and the curl operation acting on the rotational (in- 
compressible) components of velocity, VS"'(= v' ), do 
not commute. (Note that the velocity field also con- 
tains a potential (compressible) component, VS^, which 
is regular and not associated with vortices.) First, we 
shall show that V and (d/dt) commute when acting on 
S"'(= tan~^[{y — yo{t))/{x — xo{t))]). Let the instanta- 
neous velocity vector w(= dxo/dt) be chosen in the x- 
dircction, tq = and ^ be the angle around r = 0. The 
equality 

VidS'^/dt) = {w/r^)[cos{2(j))y - sin(20)i] = d{VS'^)/dt 

(12) 

holds near the vortex, ensuring that V and d/dt com- 
mute. Hence Eq.(2) is a valid equation derived from the 
Schroedinger equation even in the presence of rotational 
component. Second, the quantity dS^/dt = — sin((/))/r 
is a single- valued function, and therefore V{dS^ /dt) is 
curl-free. So will d\7S^/dt be curl- free. Upon taking a 
curl operation on Eq.(2), both sides are identically zero. 

On the other hand, we have the vortex equation of 
motion 



ing with the fact that V x {dWS^/dt) = 0, we have the 
commutation relation 



[Vx,5/Si]VS"' = w • V(V X VS"^), 



(14) 



dW X WS"- 
dt 



-w(t) • V(V X 



(13) 



where the vorticity V x VS"" = 6{r — ro{t)). Combin- 



which is non- vanishing. This completes our proof. 

From the viewpoint of dynamical evolution, we find 
that the fluid formulation of Quantum Mechanics, 
Eqs.(l) and (2), are inadequate when vortices exist in 
the system. A third equation must be supplemented, 
with the above commutation relation, Eq.(14), or the 
vortex equation, Eq.(13). This new equation updates 
the positions of vortices, and the instantaneous rota- 
tional velocity \7S^ outside the vortex can be determined 
by the Biot-Savart law of induction. Since the vortex 
velocity w, Eq.(9), depends also on the potential com- 
ponent of the flow VS'P, the rotational component VS^ 
is thereby coupled to the potential component. On the 
other hand, the potential component V5p is coupled to 
the rotational component through the momentum equa- 
tion, Eq.(2), since VS in this equation contains both \/S^ 
and V5^. It thus follows that these two topologically dif- 
ferent components of flow are non-trivially coupled. 



III. TIME-DEPENDENT SOLUTIONS OF 
QUANTUM TURBULENCE 

In order to understand the development of vortex cre- 
ation and the subsequent evolution in quantum turbu- 
lence, we solve the time-dependent solutions of the free- 
particle Schrodineger equation in a periodic box of 2048^ 
grids for 2D and of 1024^ grids for 3D. The solution can 
be given analytically through the Fourier-mode propaga- 
tor, and the periodic computational box used here serves 
only to illustrate the turbulence solutions. 



A. Initial Conditions and Potential Flow 
Turbulence 

We deflne the potential (compressional) flow in the 

Fourier space as Vk^ = fcfc • Vk. The rotational (incom- 
pressible) flow is then defined as v'' = v — v''. We choose 
the initial condition such that p = 1 throughout the en- 
tire domain and let the phase be randomly distributed 
only on the large scale. Specifically, the Fourier compo- 
nents of 5* obey a Gaussian random distribution centering 
at |k| = with the variance |Ak| equal to 20 x 27r/L for 
2D and 10 x 27r/L for 3D, where L is the domain size. 
Thus the initial flow energy is confined on large scales. 

We construct the time-dependent solution of the linear 
Schroedinger equation by using the propagator in Fourier 
space. Each Fourier component of mode k is propagated 
forward in time by multiplying the phase factor e^*^ 
to the initial Fourier mode. The updated wavefunction 
is then used to construct the fluid variables, p and v, for 
exploration of quantum fluid turbulence. 
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result of uncertainty principle. The Fourier component 
Avfcl ^ for 2D and |Avfe| ~ fc^^ for 3D, yielding a 
spectral density ^ k^^ for 2D and k^'^ for 3D. Note that 
the development of this power-law behavior is through 
cascades, where modal energy is transferred from large 
scales to small scales by local coupling in k space. This 
nonlinear coupling is similar to cascades in classical fluid 
turbulence. 

To gain a better picture of what is going on at the onset 
of vortex creation, we shall turn to the real-space flow. 
Figure (2) shows the local rise of 2D potential-flow ve- 
locity due to a large local gradient in quantum potential, 
corresponding to the situation where R = and / = 
surfaces are locally tangent to each other. Shortly after, 
two streams of counter-rotating flow develop about the 
two sides of the peak potential flow velocity to form a 
tightly bound vortex pair, corresponding to the critical 
crossing of the two surfaces. The two counter-rotating 
flow streams grow in strength until their circulations are 
quantized before the two vortices can separate from each 
other to become isolated vortices. In 3D, vortex creation 
is through formation of a vortex loop of infinitesimal ra- 
dius. The physical picture is similar to what happens 
for a vortex pair in 2D. Not until the rotational velocity 
grows to the strength where the circulation is quantized 
can the vortex loop expand to a finite radius. 



B. Vortex Turbulence 



FIG. 1: Evolution of power spectra for potential and rotation 
flows before the onset of vast excitation of rotational flows for 
2D (upper panel) and 3D (lower panel). Cascade of potential 
flow energy from large scales to small scales is apparent as 
the small-scale power progressively rises. In the mean time, 
increasing numbers of vortex pairs are forming during this 
stage, signified by the small scale excitation in the rotational 
flow power spectrum. 



Due to the nonlinearity of the fluid variable, the flow 
energy is able to couple from large scales to small scale. 
At some point in time, the smallest scale features in the 
computation domain are fully excited by the nonlinear 
coupling and the vortex flow suddenly appears, as shown 
in the velocity spectrum (Fig.(l)). It is noted that near 
the onset of the vortex flow, the small-scale potential 
flow is fully excited up to the smallest available scale and 
the spectrum of the potential flow develops a power law 
behavior, k~^ for 2D and k^'^ for 3D. This behavior lasts 
only for a short moment before the rotational flow begins 
to get excited. 

We shall pause here and offer an explanation for the 
power-law spectrum of the potential flow. In a quantum 
system studied here, the only dimensional constant is the 
Planck constant, which has dimension {length){velocity) . 
It then follows that the velocity difference |Av| of the 
potential flow over a distance |Ar| scales as |Ar|^-'^, as a 



Once created, vortices of opposite signs can move away 
from each other in 2D and vortex loops can expand in 3D. 
Starting from our relatively smooth initial condition, af- 
ter the first vortex pair is created, many more vortex 
pairs can soon be generated, tapping the free energy in 
the potential flow. In Figs. (3) and (4), we show a snap- 
shot of point vortices and vortex lines in our 2D and 3D 
computation domains. 

The turbulence spectrum long after the onset of vor- 
tex fluctuations is shown in Fig. (5). The flow eventually 
is populated with vortex fluctuations and the spectrum 
approaches a steady state that also largely obeys a power 
law as k~^ in both 2D and 3D cases. (The high-Zc upturn 
in the spectrum is an artifact of the finite-size grid and 
will be addressed later.) Despite the potential flow tur- 
bulence is suppressed in small scale after the rotational 
flow is vigorously excited, the potential and rotational 
flow turbulence energies are actually in equi-partition. 
Any vector fleld can be decomposed to the potential com- 
ponent and rotational component as V0 and z x Vt/", 
respectively, in 2D, but as V0 and A, constrained by 
V • A = 0, in 3D. That is, a 2D flow has two degrees 
of freedom arising from the independent (f) and ip, and a 
3D flow has 3 degrees of freedom given by one ip and and 
two independent components of A. If these components 
are non-correlated, equi-partition between potential flow 
turbulence and rotational flow turbulence is naturally ex- 
pected. In 2D, the equi-partitioned ratio of potential flow 
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FIG. 2: The real-space flow configuration during vortex-pair 
creation in 2D. Upper panel shows the sole existence of poten- 
tial flow (red) prior to the pair creation, where the potential 
flow locally develops a large velocity gradient. Lower panel 
shows the birth of a vortex pair, where counter-rotating flows 
(blue) grow in strength. Not until the rotational flow velocity 
grows to the magnitude where the circulation is quantized can 
the pair separate as isolated vortices. 



energy to rotational flow energy is 1, and in 3D the ra- 
tio becomes 1/2. These behaviors are what have been 
observed in our 2D and 3D cases. 

Return to the turbulence spectrum. We note that vor- 
tices are line objects when viewed in 3D. The fluid ve- 
locity is dominated by local 2D flows around the vor- 
tex lines. Therefore when the vortex lines are sparsely 
distributed in space, the vortex flow velocity has strong 
short-distance correlation transverse to the line but long- 
distance correlation along the line, leading to similar 2D 



FIG. 3: Real-space image of point vortices in steady-state 2D 
turbulence. These point vortices are captured by plotting the 
inverse density of extremely large value. The color code is in 
logarithmic scale. 



FIG. 4: Real-space image of tangled vortex tubes in steady- 
state 3D turbulence. These vortex tubes are captured by 
plotting the constant-density surfaces of very low density. 



and 3D velocity spectra in high k. In fact this spec- 
trum may be derived from the spectrum of a single vor- 
tex. The transverse velocity around a vortex |vj_| scales 
as |r|~^. Hence its Fourier component scales as k~^, and 
the spectrum as k~^. In reality, there are a number 
of point vortices and vortex loops in 2D and 3D. (See 
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FIG. 5: Steady-state power spectra of 2D (upper panel) and 
3D (lower panel) turbulence. The total flow, potential flow 
and rotational flow spectra are plotted separately. The origi- 
nal spectra have unphysical upturns at the highest wavenum- 
ber due to the grid-size singular flow velocity of vortex that 
cannot be faithfully captured by the grids. Upon removal of 
these unphysical pixels, the spectrum shows the almost per- 
fect power law in both 2D and 3D turbulence. They also 
exhibit energy equi-partition on small scales. 



Figs. (3) and (4) for real-space 2D and 3D images of vor- 
tices.) If these vortices are weakly correlated, the total 
spectrum is simply the sum of all individual vortex spec- 
tra and hence will still obey the power law. Indeed, 
since the 2D point vortex occupies a 2D area of mea- 
sure zero, so does the 3D loop vortex a zero volume, the 
vortex- vortex two-point correlation function is therefore 
overwhelmingly dominated by the zero-separation contri- 
bution, and the resulting spectrum remains to be k^^ . 

Nevertheless, the finite-separation contribution to the 
vortex-vortex two-point correlation should still exist in 
this system. It is simply due to the strong zero-separation 
correlation that the finite-separation correlation becomes 
not obvious. For example, the 2D vortex pair should have 
kept some memory of its counterpart after pair separa- 
tion, and therefore they should be correlated. In 3D the 



vortex-loop element should also be correlated with other 
elements in the same loop. Correlation of different loops, 
once exists, should be an indication of non- trivial dynam- 
ics in action. 

To uncover this weak finite-separation correlation, we 
treat every point vortex in 2D as either a signed point or 
an unsigned point and examine the two-point correlation 
function of these points. It is found that for the un- 
signed case, the vortex correlation ^(Ar) does not show 
any power-law relation with Ar, where Ar is the sepa- 
ration of two points. For the signed case, the two-point 
correlation r7(Ar) does however exhibit the screening ef- 
fect, where at small separation the correlation is negative 
and finite but the correlation abruptly drops to zero at 
some critical separation. In fact the signed vortex corre- 
lation can be fit by a Gaussian. The critical separation, 
or the variance of the Gaussian, represents the mean vor- 
tex separation, which is determined by the initial flow en- 
ergy with larger flow energy exciting more vortices and 
therefore yielding smaller vortex mean separation. Fig- 
ure (6) shows such 2D two-point correlation of signed and 
unsigned cases. 

For 3D finite-separation correlation , we also consider 
the vortex line as either a directional line or a non- 
directional line. The directional line correlation is de- 
fined as 77(Ar) = (c?li(r) ■ c?l2(r + Ar)), and the nondirec- 
tional line correlation as i^(Ar) = (d^i(r)(i/2(r + Ar)). It 
turns out that the nondircctional line-line correlation ^ 
obeys a universal Ar~^ power law. Although the direc- 
tional line-line correlation rj also obeys Ar~^ power law 
at small separation, it abruptly drops to zero at a finite 
separation. The directional line-line correlation rj reveals 
the average size of the vortex loop at the separation when 
the correlation drops to zero. Figure (7) depicts the 3D 
line-line correlation. Comparing these correlation func- 
tions in 2D and 3D cases, we find that the dimensionality 
plays a role in the final states. Two-dimensional quantum 
turbulence in our model appears farther from a relaxed 
state than the 3D case since small-scale mixing in 3D 
turbulence is well in place in the final steady state that 
yields a power-law correlation. 

We now turn to the upturn at the highest-fc in the 
energy spectrum. It is actually caused by numerical er- 
rors in the finest grid at the vortex. The velocity field is 
singular at the vortex, whose exact location can be any- 
where within the finest grid, giving unequal weights to 
the grids immediately surrounding the vortex. When the 
vortex is much closer to a particular grid, the velocity 
at that grid will be much larger than other nearby grids. 
Therefore, this grid will appear as an erroneous Dirac 5- 
function, thereby giving an upturn with the first power of 
k in the spectrum at the smallest scale. We suppress this 
problem by limiting the magnitude of velocity to an ap- 
propriate upper bound so as to eliminate the contribution 
of the said (5-functions. (Only 20-30 grids are processed 
among a total of 2048^ grids that contains about 100 
point vortices in the 2D case, and the 3D case also has 
a similar ratio of processed to unprocessed grids.) Fig- 
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FIG. 6; Two-point correlation functions of signed point vor- 
tices rj and unsigned point vortices ^. The signed vortex cor- 
relation rj can be fit by a Gaussian and the variance of the 
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FIG. 7: Line-Line correlation of directional vortex lines rj and 
non-directional cortex lines ^. Both follow a power law Ar~^ 
at small separation. 



ure (5) also shows the spectra after this regularization 
procedure. We emphasize that such a grid-scale error 
is unavoidable since the local physics is singular; such 
an error can accumulate over time when evolving a non- 
linear schroedinger equation numerically. Only in this 
soluble model adopted here can the error in the solution 
be controllable and does not accumulate. 

Quantum vortices are long-lived coherent structures in 
dissipationless turbulence. Once created, they can only 
come to destruction through pair annihilation in 2D and 
through shrinkage of ring vortex in 3D. When vortex 
lines are sparsely distributed in space, vortex collisions 
are rare events, thereby making them long-lived objects. 



Vortex-line reconnection arises in 3D from vortex col- 
lision and is a more common phenomenon than vortex 
annihilation. Details of 3D vortex-line reconnection are 
depicted in Fig. (8). Also clear in Fig. (8) is the existence 
of helical interference fringes along the cores of 3D vortex 
lines. These fringes appear to vary rather slowly and are 
the TO (mode number around the tube)= 2 density fluc- 
tuations of varying n(mode number along the tube). We 
mentioned previously that the density profile cut across a 
vortex-tube is generally elliptical near the null, and this 
is consistent with the m = 2 density fluctuation observed 
at the 3D vortex core. We also note that Kelvin waves are 
the TO = 1 fluctuations and oscillate due to the " tension" 
of a vortex line. Thus, Kelvin waves appear qualitatively 
different from these to = 2 fringes. 

In 2D quantum turbulence, inverse cascade has not 
been observed in our system. We explore this issue by 
having the energy of initial potential flow concentrating 
on the mid scale, hoping to provide an ample space on 
the large scale for inverse cascades to take place. How- 
ever, the large-scale flow energy in fact does not change 
throughout the evolution, indicating the non-existence of 
inverse cascade in this model. Treating unsigned vortices 
as individual points, we again examine the vortex-vortex 
correlation, aiming to find out whether these points have 
any tendency to cluster when these point vortices are 
created on an intermediate scale. Like the situation ex- 
plored previously, we find that 2D vortices do not develop 
long-range correlation. This aspect is in great contrast to 
classical 2D fluid turbulence, where vortices of the same 
sign tend to aggregate to form large vortex clusters. 



IV. CONCLUSIONS 

By exploring the exactly soluble Schroedinger equation 
for free-particles, we have derived the following generic 
results in this work. These results can also, or has been 
observed to, arise from a nonlinear Schroedinger system. 
(I) Quantum vortex has a axisymmetric rotational ve- 
locity but a non- axisymmetric potential velocity, both of 
which diverge near the vortex as r~^. (2) Whenever vor- 
tices appear in the system, the density and momentum 
equations of the fluid representation will become inade- 
quate to evolve the system, and it requires a third equa- 
tion, the vortex equation of motion, for the full descrip- 
tion of evolution. (3) Unlike classical vortex, quantum 
vortex is a wave that is not frozen into the matter and 
it propagates at a velocity different from the local mat- 
ter velocity. (4) The motion of quantum vortex is not 
entirely induced by other vortices and hence its propaga- 
tion velocity is different from the rotational component 
of local matter velocity. (5) Equip-partition between ro- 
tational flow energy and potential flow energy can be 
established on small scales in fully developed quantum 
turbulence. 

Specific results for our 2D and 3D free-particle mod- 
els are that (a) starting with an initial potential flow. 
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FIG. 8: 3D reconnection of vortex loop, which is the major 
activity for the maintenance of steady-state quantum turbu- 
lence. The evolution proceeds from the upper to bottom pan- 
els. We also note the small-scale helical fringes in existence 
along the cores of vortex lines. These helical fringes appear at 
the very core of a vortex tube, as they can only be captured 
when the constant-density surfaces assume an extremely low 
value, lower than that used to plot Fig.4. 



turbulence cascades toward small-scale is observed prior 
to rotational flovifs are excited, and near the transition 
to the excitation of rotational flow the energy spectra 
of potential flow turbulence in 2D and 3D obey differ- 
ent power laws; (b) the energy spectra of fully developed 
2D and 3D turbulence with vortex excitation both obey 
the same power law, which corresponds to the power 
spectrum of a single vortex; (c) Ar~^ power-law vortex- 
vortex correlation is observed only in 3D but not in 2D 
turbulence; (d) there is no evidence of inverse cascade in 
2D turbulence; (e) there is no evidence of Kolmogorov 
energy cascade in either 2D or 3D turbulence. 

By comparing our general and specific results, it is ev- 
ident that small-scale features derived from this linear 
Schroedinger equation can be in common with those de- 
rived from a nonlinear Schroedinger equation. But the 
large-scale features, such as vortex-vortex correlation, 
cannot, because the wave nonlinearity and dissipation, 
which are lacking in this linear model, must play an es- 
sential role in shaping the global dynamics. Along this 
line, we note that the steady-state turbulence spectrum 
is maintained not by the conventional cascade processes 
in this model, since no dissipation is involved. In fact, en- 
ergy cascade is only observed in the initial potential flow 
phase as a transient before vortices are excited. Once 
vortices are fully excited, the steady-state energy 
spectrum found in this model only reveals the spectrum 
of a collection of individual vortices and contains no in- 
formation about correlation among different vortices. 

Whether it is the dissipation or the wave nonlinear- 
ity that is the prime factor for producing the long-range 
vortex- vortex correlation is an important issue for under- 
standing the actual quantum turbulence. In fact, in the 
presence of small-scale dissipation the steady-state 
spectrum must fail, as it would have yielded a diverging 
dissipation rate. Whether the Kolmogorov spectrum may 
be established in a dissipative linear Schroedinger system 
remains to be investigated. In the case of a nonlinear 
Schroedinger system, Kelvin wave cascades are found to 
play a central role in coupling to dissipation [8,9]. On 
this particular issue, our observation for a vortex loop of 
a linear Schroedinger system is that it can undergo large- 
amplitude loop distortion so that the loop is self pinched 
and reconnects to form two smaller loops (c.f.. Fig. (8)). 
This reconnection process can couple to dissipation effi- 
ciently if some small scale dissipation is included in the 
system. On the other hand, despite Kevin waves also ride 
on a vortex loop and can give rise to large-amplitude loop 
distortion, the picture of Kelvin wave cascade from long 
waves to short waves and finally to dissipation on a vor- 
tex loop is quite different from our observation of vortex 
loop dynamics of a linear Schroedinger system. 

Finally, we should mention that the present dissipa- 
tionless linear system in a periodic box can be recurrent 
in a finite time when the longest wave in the system oscil- 
lates over one period. The steady state turbulence pre- 
sented in this work is found to be already established 
within a small fraction (10%) of the recurrent time. The 
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Poincare recurrence time of this model is proportional 
to the square of the periodic box size, L^. Despite the 
dynamics of this system may appear chaotic, it is not 
quite a chaotic system due to the algebraic, instead of an 
exponential, dependence of the recurrence time on the 
degrees of freedom. 
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